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Measurements of the persistent current in a ring containing a quantum dot would afford a unique 
opportunity to finally detect the elusive Kondo screening cloud. We present the first large-scale 
numerical results on this controversial subject using exact diagonalization and density matrix renor- 
malization group (RG). These extremely challenging numerical calculations confirm RG arguments 
for weak to strong coupling crossover with varying ring length and give results on the universal 
scaling functions. We also study, analytically and numerically, the important and surprising effects 
of particle-hole symmetry breaking. 
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The screening of an impurity spin by conduction elec- 
trons, the Kondo effect, is believed by many to be associ- 
ated with the formation of a "screening cloud" around the 
impurity with a size £k = vf/Tk where vf is the Fermi 
velocity and Tk is the Kondo temperature, the charac- 
teristic energy scale associated with this screening 0. 
This fundamental length scale, which has been called the 
"holy grail" of Kondo research |2J, and which from the 
above estimate can be as large as 1 micron in typical sit- 
uations, has never been observed experimentally and has 
sometimes been questioned theoretically. Although tra- 
ditionally associated with dilute impurity spins in met- 
als, the Kondo effect has been observed more recently in 
nano-structuresiSii. The electron number on semi- 
conductor quantum dots, weakly coupled to leads, can 
be varied in single steps with a gate voltage. When this 
electron number is odd, the quantum dot generally has a 
spin of 1 /2 and can act as a Kondo impurity, screened by 
electrons in the leads. By attaching the dot to quantum 
wires of length L, the Kondo screening cloud size could be 
measured from the dependence of conductance properties 
on L 00] A particularly simple case is when the dot is 
attached to a single quantum wire which forms a closed 
ring. Two simple tight-binding models have been consid- 
ered corresponding to an "embedded" or "side-coupled" 
quantum dot (EQD or SCQD). Suppressing electron spin 
indices, the corresponding Hamiltonians are: 
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where H K = JkS ■ (c{ + c} / _ 1 )|(ci + cl-i), and 
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where Hk = JkS -Cq^co, and Cj„ annihilates an electron 
at site j of spin er; S a are S=l/2 spin operators. We 



generally set t = 1 in what follows. A magnetic flux is 
added to the model by adding appropriate phases to the 
hopping terms and the resulting persistent current, j, is 
measured. 

The EQD must have j = when the Kondo coupling, 
Jk = 0, since then the sites 1 and (-1) are decoupled 
from each other. On the other hand, the SCQD has a 
large j when Jk = since the Hamiltonian then reduces 
to that of a free ring with periodic boundary conditions. 
When the Kondo coupling becomes large (Jk ^> t) es- 
sentially the inverse behavior occurs, due to the forma- 
tion of a Kondo screening cloud. In the strong coupling 
limit of the EQD the screening electron goes into the 
symmetric orbital on sites 1 and -1. This allows resonant 
transmission through the anti-symmetric orbital on these 
sites and the ideal saw-tooth like j of a free ring occurs, 
when the system is at 1/2-filling. On the other hand, 
for the SCQD, the screening electron sits at site-0 and 
completely blocks all current flow in the strong coupling 
limit. The more interesting question is the behavior of j 
in these models for small Jn/t at large length, L. The 
Kondo length scale grows exponentially as Jx/t — > 0. 
The persistent current was predicted 0] to be given 
by universal scaling functions of £,k/L and the dimen- 
sionless magnetic flux through the ring, a = e$/c: 



jL/(ev F ) = /(£ic/L,a), L,£ K ^a, 



(3) 



with a the lattice constant. The crossover functions, /, 
are very different for the EQD and SCQD and also de- 
pend on the parity of the electron number in the ring, 
N. However, they are otherwise expected to be uni- 
versal in the small Jk limit, not depending on details 
of the dispersion relation, electron density, the range of 
the Kondo interaction, etc. In previous work by one of 
us and P. Simon 0, it was argued that, in the limit 
L 3> £k> strong coupling behavior occurs, since the effec- 
tive Kondo coupling is expected to become large at large 
length scales. Thus it was predicted that J^k/L, a) ap- 



2 



proaches the value for an ideal ring at £x / L — > for the 
EQD but approaches zero for the SCQD, for either par- 
ity of N . The former prediction is in disagreement with 
other studies where jL for the EQD was predicted 

to have very different L-dependence, whereas the SCQD 
prediction is in stark contradiction with the conclusions 
of other authors 0, 0] who predicted that j attains 
the value for an ideal ring, at £k/L — > 0. However, the 
results of Refs 0, are in agreement with the theoreti- 
cal, numerical (u p t o L = 8) and mean field results of 
Refs. ITR Ei l Ha I20L Other theoretical studies include 
Refs. ll3l ll4L Il5l lltl The behavior of the persistent cur- 
rent, j in quantum dot systems is therefore the subject 
of considerable controversy. 

In this letter we attempt to settle this controversy 
using exact diagonalization (ED) and Density Matrix 
Renormalization Group (DMRG). The calculations have 
been performed using fully parallelized programs on dis- 
tributed SHARCNET facilities. The DMRG work turns 
out to be very cpu-intensive due to the necessity of us- 
ing periodic boundary conditions to obtain a persistent 
current, the complex form of the Hamiltonian at non- 
zero flux, the necessity of calculating a derivative of the 
groundstate energy to obtain the current: j = —edE/da, 
and the fact that the current is o(l/L). In the case of 
the EQD we find good agreement with the expected scal- 
ing picture and obtain useful results on the crossover 
functions, /. For the SCQD we find an extremely slow 
crossover with varying L or Jk but we present both 
numerical and analytical evidence that j scales to zero 
at small £k / L at 1/2- filling, as predicted by the previ- 
ous RG approach. We show that particle-hole symme- 
try breaking leads, for small Jk/t, to small non-universal 
corrections to jL. Surprisingly, these produce a small 
non-zero value of j e L, the current for even N , at L — > 00 
for the SCQD. 

Embedded Quantum Dot: In Fig. JQ) we plot j vs. a/ir 
for a fixed Jk = 1, and various L = N. Here both L and 
N includes contribution from the impurity site. Note 
that a crossover is seen between the weak coupling be- 
havior at smaller lengths (small and sinusoidal) to strong 
coupling ideal ring behavior at the largest lengths (larger 
and the saw-tooth). To further illustrate this crossover, 
we focus on one value of the flux, a/ir = 0.05. In 
Fig. we plot jL/evp vs. L for this fixed value of 
the flux for L = N both even or odd. Small numeri- 
cal errors in the DMRG calculations with m=(1024,512) 
are visible beyond L=24 when compared to results with 
m=(1024,2048). As can be seen from Figs. (JTJ) and (0) 
the current depends strongly on whether L = N is even 
or odd. For L — N even there is a difference between 
N = 4p and N = 4p + 2, visible in Fig (J2J), that can be 
absorbed into a re-definition || of the flux a — a + irN/2. 
From the results shown in Figs. and (0 we conclude 
that jL increase with L towards the free ring limit con- 
tradicting Refs I9L Hol 
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FIG. 1: The current at 1/2-filling for N = 4p - 1 (A) and 
N = 4p (B) for the EQD. Results are shown for a number of 
system sizes with Jk = 1 as a function of a/ir. The results 
for the small system sizes are obtained using exact diagonal- 
ization methods while the larger system sizes (circles) have 
been obtained using DMRG techniques. 




L (With impurity) 

FIG. 2: The current jL/evF for Jk ~ 1 and a/n = 0.05 as 
a function of L for the EQD. The results for the small sys- 
tem sizes are obtained using exact diagonalization methods 
(•) while the larges system sizes have been obtained using 
DMRG techniques with (mL = 1024,m fl =512) (o). The □ in- 
dicate DMRG results obtained using (mi=1024, m fl =2048). 
The dashed lines indicate the free ring result. 



Weak and strong coupling results @, and Fig. ©indi- 
cate that j(a) has period 27r for N = L even but period 
7r for N = L odd. The latter result follows rigorously 
at 1/2-filling from particle-hole (P-H) symmetry which 
takes a — > 7r — a for N = L odd together with time re- 
versal which takes a — > —a. Away from 1/2-filling, when 
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FIG. 3: The current at 1/2-filling jL/ev F , at a/% = 0.4 
(TV = 4p) o/tt = 0.8 (TV odd), plotted vs. L/£x obtained by 
combining our results for various values of Jk for the EQD. 
We fix ^k(Jk = 0.3) = 1 and obtain the remaining £k rel- 
atively to £,k(Jk = 0.3) by rescaling. The inset shows the 
obtained £k{Jk) as a function of 1/Jk- 



P-H symmetry is broken, we might expect that the pe- 
riod of j(a) would be enlarged to 2ir. This can be checked 
at both weak and strong coupling. In the weak coupling 
limit, in addition to the terms calculated previously jjj, 
we find an additional term, of period 2ir in the current 
for TV odd, only present away from 1/2-filling (N/L ^ 1): 

djo ra MM sin & sin 2 ( nN /2L) 271 , ( 4 ) 

where a = a + tt(N — l)/2. In the weak coupling limit, 
the period 27r term in j for odd TV can be understood, in 
the continuum limit formulation with a linearized disper- 
sion relation and a wave- vector cut-off which is symmet- 
ric around the Fermi surface, as arising from particle-hole 
symmetry breaking potential scattering terms in the ef- 
fective Hamiltonian, of 0(J K ). These terms are strictly 
marginal under renormalization group transformations 
and hence their effects do not grow larger with growing 
L/£k, but remain of 0(J K ), where Jk is the bare cou- 
pling, at all length scales. Thus the period 27T term in j 
becomes negligible at all length scales, when L>^, in 
the limit of small bare coupling and hence do not appear 
in the universal scaling functions, f(£,x/L,a). 

We have also attempted to numerically calculate an 
approximation to the scaling function / by rescaling ED 
results out to L = 13 for a range of Jk- Our results are 
shown in Fig. where jL/evp is plotted vs. L/£k for 
a fixed a/ir = 0.40(0.80) for L = N odd(even). If 
is fixed at a given J K , £{J K ) can be estimated by rescal- 
ing L/£k until jL/evF, calculated with J K superposes 
jL/evF calculated with J K and so forth. The reults of 
Fig. @ nicely confirms the scaling picture and shows 
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FIG. 4: (A) The current jL/e at 1/2-filling for the SCQD 
for a/n = 0.15 and Jk = 0.75, 1.60, 4.00 as a function of 
L. Exact diagonalization results (solid symbols) are shown 
for small system sizes and compared to DMRG results with 
(m L =1024,m fl =2048) for L > 13 (shaded symbols). For com- 
parison a function with pure 1/L behavior is shown. Note that 
the current goes to zero in the strong coupling limit. (B) The 
current jL/e at 1/4-filling for the SCQD for a/ir = 0.75 and 
Jk = 4.00 as a function of L. Exact diagonalization results 
are shown for system even N (o) up to N = 10(L = 18) and 
for odd TV (•) up to TV = 11(L = 20). The dashed line is a 
power-law fit to the points with TV = 5 — 11, of slope=-1.00. 



that / is an increasing function of L/^k- 

Side Coupled Quantum Dot: [Now L does not include 
the impurity site while TV does include the impurity elec- 
tron. For N odd a = a + ir(N — l)/2 and for N even 
a = a + irN/2.] In Fig. (0J\) we show exact diagonal- 
ization and DMRG results for the SCQD at a = 0.15-7T 
analogous to Fig J5J. These results are consistent with 
jL going to zero at I > (r- however it appears neces- 
sary to go to extremely large L/^k to see this behavior. 
A calculation of the scaling function, /, equivalent to Fig. 
@, shows that / for the SCQD is a decreasing function 
of L/£k in contradiction to the results of Refs. fill H^. 

In order to test our conjectured scaling behavior, we 
have studied both analytically and numerically the strong 
coupling limit, Jk ^ 1 in greater detail. In this limit 
one electron sits at site and forms a singlet with the 
impurity spin, effectively cutting the ring at the origin. 
Perturbation theory in 1/Jk generates terms in a low 
energy effective Hamiltonian which couple the two sides 
of the quantum dot and thus allow for a small non-zero 
current. By doing perturbation theory to 3 1C ^ order in 
1/Jk, we obtain an effective Hamiltonian valid at Jk 3> 
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where the tunneling terms are: 
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Here rii is the electron number on site i. We have ignored 
additional terms in the effective Hamiltonian which do 
not contribute to the current at low orders in 1/ Jk- We 
may now evaluate the current, up to 0(1/ J K ), by calcu- 
lating the a-dependence of the groundstate energy in first 
order perturbation theory in Ht- Considering 1/2-filling, 
this gives, for odd or even N, 
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The absence of period 2tt terms for N even can be un- 
derstood from the P-H symmetry of the unperturbed 
Hamiltonian in Eq. (JSJ which is broken, for N even, 
by Ht- Note that these formulas predict that j goes to 
zero as 1/L 2 at large L. In the strong coupling limit, 
this amounts to an analytic proof that the current is not 
the same as for the ideal ring, in contradiction with the 
claims of Ref. [llj. We have verified that Eqs. Q are in 
excellent agreement with our numerical results for large 
Jk- However, we find that surprisingly large values of 
Jk, of about 100, are necessary before these formulas 
become accurate. This is consistent with the results of 
our study of length dependence of j which suggests a very 
slow cross over from weak to strong coupling behavior as 
L/£k is varied. 

From Eqs. Q it is also possible to calculate the current 
away from 1/2-filling, at large Jk we find: 
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Note that j L/e is 0(1/ L) and actually gets smaller as we 
move away from 1/2-filling. We sec that j Q L approaches 



zero as L — > oo. Surprisingly, we see that j e L attains a 
non-zero limit as L — > oo, away from 1/2-filling, in stark 
contrast to j L and the behavior of both j a L and j e L 
at 1/2-filling. We have numerically verified Eqs. JBJ in 
detail finding excellent agreement for Jk > 100. For an 
intermediate coupling of Jk = 4 we show ED results for 
jL/e in Fig |QJ3), at a flux a/Tr = 0.75 and 1/4-filling, 
clearly approaching a non-zero limit for N even with op- 
posite sign and a much smaller amplitude than that of 
the ideal ring. For odd N, jL/e ~ const/L in accordance 
with Eq. JSJ. This 0(1/ ' L) behavior of j e away from 1/2- 
filling results from non-universal particle-hole symmetry 
breaking terms in the low energy effective Hamiltonian. 
For small Jk/D (where D w t is the bandwidth) these 
terms are small, and remain so under renormalization. 
They contribute to j e at L ^> which has the form: 
j e = (ev F /L)[(A£ K /L)sm2a + B(J K / D) 2 sina] where 
A is a universal constant and B is a non-universal con- 
stant, both of 0(1). Thus the universal (1/L 2 ) behavior 
is destroyed at large length scales, « £,k(L> / ' Jk) 2 ^ £k, 
away from 1/2-filling, in contrast to the claim in Ref. 
(@). Note however that the current at these large length 
scales is smaller by a factor of (Jk/D) 2 than that of an 
ideal ring, in contrast to the claim in Ref. HT 



Following NozieresQ we have developed a local Fermi 
liquid theory description of the low temperature fixed 
point of the SCQD. This fixed point corresponds to the 
impurity forming a singlet with one conduction electron 
and the remaining low energy electrons being repelled 
from the origin, corresponding to a broken ring with 
zero persistent current. In the 1/2-filled case, the lead- 
ing irrelevant operators at this fixed point, which con- 
trol the persistent current at large L, are (ip+aip-) 2 , 
ipXaif}- ■ (tp' ! + aip+ + ip_aip-) and their complex conju- 
gates. Here + and — label the electron operators on the 
right and left hand side of the impurity. In the limit of 
weak bare Kondo coupling, or large the coupling con- 
stants in front of these terms can all be fixed uniquely up 
to one overall factor with dimensions of inverse energy, 
proportional to the inverse of the Kondo temperature. 
This theory predicts that the persistent current scales to 
zero as eVp/TKL 2 at 1/2-filling, or for even N at arbi- 
trary filling, as well as making various other predictions 
that could be compared with numerical simulations and 
experiments. |2lj 



In conclusion, we have given strong numerical evidence 
and analytical results in support of the scaling picture of 
the persistent current in quantum dot systems resolving 
a number of outstanding controversies. This research is 
supported by NSERC of Canada, CFI, SHARCNET and 
CIAR. IA acknowledges interesting conversations with P. 
Simon. 
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